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ABSTRACT 


LIMIT  DISTRIBUTIONS  OF  KOLMDGORDV-SMIRNOV 
TYPE  STATISTICS  UNDER  A FIXED  ALTERNATIVE 
WITH  ESTIMATED  LOCATION  AND  SCALE  PARAMETERS 

Much  attention  has  baan  devoted  to  Monte  Carlo  simulations  of  the 
power  of  Kolmogorov-Smirnov  type  goodneee-of-f it  statistics  when  nuisance 
parameters  of  the  hypothesised  distribution  sre  estimated.  Here  we  con- 
sider the  asymptotic  behavior  of  such  statistics  at  fixed  alternatives 
when  location  and  scale  parameters  are  estimated.  It  Is  shown  that  suit- 
ably normalised  Kolmogorov- Smirnov  statistics  converge  In  distribution 
to  Gausslan-related  random  variables  depending  on  the  alternative  distri- 
bution and  the  maximum  deviation  between  the  null  and  alternative  dis- 
tribution functions.  The  work  of  Raghavachari  (1973)  la  thus  extended 
from  simple  hypotheses  to  the  case  of  composite  hypotheses  with  estimated 
nuisance  parameters. 
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^ ' Introduction . Lot  X^,  ■ • l(  bo  o sample  of  lndopondont  obser— 
votiona  on  a continuous  distribution  function  P.  Consider  the  oompoeits 
null  hypothesis 


V •>»  for  •m*  8 e 0, 

whore  C is  s specified  cdf  and  6 is  a specified  parameter  space  for  the 
nuisance  parameter  6.  A test  of  Hq  may  readily  be  constructed  from  the 
sample  distribution  function  and  a consistent  estimator  0 of  6.  In  partl- 
cular,  with 

Un<t;  G»  V “ n"1J1IlG(XiJ  V * tl,  0 S t s 1, 

idiere  I IE]  denotes  the  indicator  of  the  event  B,  a measure  of  discrepancy 
between  the  observations  and  the  hypothesised  distribution  is  given  by 
wodified  empirical  stochastic  process 

Vn(t;  G)  - n*[Un(t;  G,  0Q)  - t],  0 £ t £ 1. 

Many  test  statistics,  including  the  Kolmogorov- Smirnov  statistics  with 
estimated  parameters,  can  be  written  as  functionals  of  the  stochastic 

process  V^C* ; G).  Specifically,  the  one-sided  Kolmogorov- Smirnov  statis- 
tics are 

Dn+  ‘ V'1  *>•  V'1"*  ».<*!*>. 

ostsi  Ostsi 

and  the  two-sided  K-S  statistic  is 

D - sup  IV  (t;  G)|. 

Ostsl  n 


* 
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Stephens  (1974)  gives  a comprehensive  treatment  of  statistics  which  can 
be  given  similar  representations.  A more  extensive  class  of  test  sta- 
tistics can  be  shown  to  be  asymptotically  equivalent  to  such  functionals. 

See  Green  and  Hegazy  (1976)  end  Kumar  and  Pathak  (1977). 

The  weak  convergence  properties  of  the  modified  empirical  process 
and  corresponding  test  statistics  have  been  extensively  Investigated  both 
under  the  null  hypothesis  and  under  sequences  of  oontiguoua  alternatives 
by  Durbin  (1973),  Neuhaus  (1976),  and  Hood  (1978).  However,  power  stu- 
dies of  statistics  based  on  V (•;  G)  have  beam  restricted  mainly  to  Monte 

n 

Carlo  simulations. 

Here  we  consider  the  asymptotic  behavior  of  Dq+  d~  R 1fn^Tr 

* n n 

a fixed  alternative.  In  particular,  restricting  attention  to  location 
and  scale  families  of  distributions  under  both  the  null  and  the  alterna- 
tive hypotheses,  we  extend  the  results  of  Raghavachari  (1973),  who 
treated  the  case  of  simple  fixed  alternatives,  to  the  case  of  estimated 
nuisance  parameters.  Section  2 presents  the  main  results.  Section  3 the 
proofs,  and  Section  4 an  illustration  testing  normal  versus  Cauchy. 

2.  Results.  Consider  the  goodnes3-of-fit  problem  of  testing 

Hq:  F(x)  - G((x  - aG)/3G>»  -«  < x < «,  for  some  (Og,  8g),  0g  > 0, 
versus 

Hj_s  F(x)  - H((x  - a^/Bg),  -»  < x < «,  for  some  (e^,  BH),  0H  > 0. 

The  approach  In  Section  1 is  motivated  by  the  well-known  result  that  under 
Hq  the  random  variables  G(X^;  (Og,  0g))»  1 * 1 s n,  are  Independent  uni- 
form (0,  1)  variates.  In  this  case  the  sample  distribution  function 
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_»  n 

un(t;  G*  <<V  " n I IJc(*ii  (oqi  6q))  stl.Ottil, 

la  strongly  conalstant  astlnator  of  tha  uniform  distribution  function 
U(t)  • t,  0 i t s 1.  Tha  K-S  statistics  can  than  ba  vrlttan  «•  func- 
tionals of  tha  anplrlcal  stochastic  process 

«n(t;  C)  - n*lUn(t;  C,  (og,  Bc))  - t],  0 * t s 1. 

That  la,  sup  V (t;  G),  Inf  W (t;  G),  and  sup  |W  (t;  G) | . 

0«*1  n Oitsl  n Ostsl  “ 

Than  weak  convergence  of  Wq(.;  g)  to  tha  tied -down  Hlanar  process  W° 

yields  tha  Halt  distributions  of  tha  K-S  statistics.  (See  Billingsley 

(1968)  for  details.)  Hsre  W°  la  tha  Gaussian  process  determined  by 

8{V°(t) ) - 0.  0 s t i 1,  and  B{W°(s)W°(t))  - nln(s,  t)  - st,  0 S s,  t S 1. 

Similarly,  to  test  Hq  as  given  above,  we  replace  Wn(»;  G)  by  Its 

analogue  V (•;  G)  based  on  estimates  (a  , 8 ) of  the  nuisance  parameters 

(°G»  ®g*  * ***  r*9ulr*  thet  (on,  6^)  be  consistent  estimators  satisfying 

conditions  which  Insure  the  weak  convergence  of  Vr(* ; g)  under  H^.  See 

Aesuaptlons  A and  B(l).  These  conditions  also  Insure  that  G(X, ; (o  , $ )) 

1 n n 

converges  In  probability  to  a uniform  random  variable  and  that 
Un(t;  G,  (oQ,  §n))  is  a consistent  estimator  of  U(t),  0 < t s 1. 

Hoover,  under  H^,  (an,  Sr)  need  not  converge  to  t he  appropriate 

location  and  scale  parameters  of  H but  we  aeaume  the  existence  of  cl,  and 

* 4 . * 

BK  > 0 such  that  n ($n  - c^)  - 0p(l)  and  n (Bn  - BR)  - 0p(l) . (In  the 

sequel,  all  probability  statements  will  refer  to  the  alternative  hypo- 
thesis, unless  otherwise  Indicated.)  Then  H((X1  - 4 )/$n)  will  converge 
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in  probability  to  H((X^  - o^)/^),  which  any  not  bo  unifora.  In  ordor 
to  produco  asyaptotlcally  unlfora  variates,  tho  probability  Integral 
tranafora  X((X^  - S^)/^)  la  needed,  where 

X(x)  • HKggX  ♦ (o^  - oH))/eK],  -•<*<•. 

In  terse  of  thla  tranaforaatlon 

D+  - aup  (V  (t;  X)  ♦ n*[t  - oT^t)]}, 

° Oitsi  n 

with  D~  and  Dq  repraaantad  analogously. 

To  Insure  the  weak  convergence  of  Vq(*;  K)  under  H^,  we  as suae  that 

X satisfies  Asauaptlons  A and  (a_ , $ ) satisfies  Aasuaptions  B,  as  follows. 

n n 

A»»wf>tion»  A. 

(I)  X'  is  positive  on  the  support  of  X; 

(II)  X"  la  bounded  and  xK'(x)  ♦ 0 as  |x|  ♦ 0.  0 

Aimption*  P. 

(1)  There  exist  randoa  variables  Y and  Z such  that  for  each 

k • 1,  2,  ...  and  0 < tj tfc  < 1,  the  vector  W°(tj),  ...,  W°(tk) , Y,  Z) 

la  Multivariate  noraal  with  seen  vector  0 and  (k+2)*(k+2)  covariance 
aatrix  S . and  la  the  weak  convergence  Halt  of  the  randoa  vector 

0,n^ti»  Wt  •*..  VV  10 » “*(®n  * V^X*  n%(®n  * *X>/$X>; 


(11)  There  exist  conatanta  c^,  c^ , 6^  and  with  aln(6^,  4^)  > 0 
such  that  for  all  n sufficiently  large 
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Jr  . -(!♦«.) 

- c^|/8g  > A}  < CjA 

and 

k . -U+6J 

P{n*l#n  ’ * X)  K C2X  • D 

With  these  uiuftlooi,  Vr(*  ; K)  converges  weakly  to  the  Gaussian 
stochastic  process 

V°(t)  - W°(t)  + K'flf^tW  + K’(K“X(t))K_1(t)Z,  0 $ t n, 

(see  Wood  (1978),  (2.2)),  while  n^[t  - GK~X(t)]  Is  unbounded  for  G * K 

L .1 

and  Is  dominated  by  n^aup  It  - GK  (t)|.  Setting 

0*tsi 

•+  - sup  It  - GK“1(t)l,  8"  - inf  (t  - Gk’l(t)l 
OSt*l  Ostsi 

sad 

©*  • (t * t - GK-1(t)  - •*■},  e"  - (t:t  - GK-1(t)  - e"), 

if  we  can  show  that 

»:  - ««p  Jv  <t;  *>  ♦ «*It  - Gr"1(t)))  ♦ 0(1) 
n tee*  “ P 

(2.1)  • sup  .V  (t|  K)  ♦ + o (1) 

t«e*  B P 

and  slallarly 

(2.2)  D‘  - Inf  V (t;  K)  ♦ n*8~  ♦ 0(1). 

“ t«©~  “ p 

we  will  have  froa  the  alaoat  sure  continuity  of  the  saaple  paths  of  V° 
that 


TUORBM  1.  Under  Aaewyticne  A end  B,  for  every  o. 


and 


li>»,P{D*  - nV*  so}  - P{aup  . V°(t)  s o) 

t*e 

llnP{D~  - nV  so)  - P{lnf  V°(t)S  o ). 

tee 

Typically,  0+  and  6 consist  of  on#  point  sach,  In  which  cass  tha 
Halt  distributions  ara  normal. 

Similarly,  setting  6 - ■a«{e'f,  -8~},  Q1  - (t«  t - Gr”1(’ ) - • ), 

®2  " t - GK  *(t)  - -9},  and  9 - u ©2,  If  wa  can  show  that 

(2.3)  Dn  - •uP^glVjjtt;  K)|  ♦ n*8  ♦ op(l)  , 

wa  will  hava 

THEOREM  2.  Under  Aeewnptione  A and  Bt  for  soarty  a, 

14»^,P(D  - n*0  so}  - P{  max  [sup  V°(t),  - Inf  V°(t)]  S a). 

tcOj  t«e2 

(Hots  that  ona  of  8^  or  8^  nay  ba  anpty.  In  such  a casa,  wa  «dopt 
tha  convantlon  that  tha  sup  over  tha  anpty  sat  Is  -•».) 

3.  Proofs.  In  ordar  to  coaplata  tha  proof  of  Thao ran  1,  wa  naad 
only  show  (2.1)  and  (2.2).  Slnca  tha  argunsnts  In  each  casa  ara  Idantlcal, 
wa  will  coosldar  only  (2.1). 

By  construction  of  K and  by  Assumptions  A and  B(l) , It  follows  froa 
Wood  (1978) , formulas  (3.3)  snd  3.1A),  that 

•up  IV  (t{  K)  - A (tj  K)l  - o(l), 

OStsi  *»  “ P 

whara 


An(ti  K)  - n%{Wn(t;K)  ♦ K,(K"1(t))I(dn  - o^)/BK  ♦ lTl(t)(Sn  - Bj)/^]). 
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for  0 s t s 1.  Thus  (2.1)  Is  equivalent  to 

aup  U (t{  K)  ♦ n*(t  - af^Ol)  - aup  .A  (t;  K)  + nV  ♦ o (1). 
Ostsl  ° tce+n  P 

Floor  OF  THEOREM  1.  For  brevity,  we  put 

B*  - aup  .{A  (t;  E)  ♦ n*(t  - af1^)]) 
n t€0+  n 

• aup  ,{A_(t;  K)  + n*6. 
t<6 

Next  for  every  Integer  lt  A - als«^i  A^},  and  0 s y s 1,  define 

(3.1)  S(y,  *)  - ( tiOSUl,  |GE_1(y)  - GE_1(t)|  < e-*5) 
and 

(3.2)  S - (S(t,  t)  s t c e*}. 

Since  la  continuous , 0+  la  conpact  and  S,  being  an  open  cover  for 
0+,  baa  a finite  cover  T(T  c $).  By  the  leans  given  at  tha  conclusion 
of  this  section,  T has  cardinality  Y not  exceeding  2[[2^6]).  where 
11*1)  denotes  greatest  integer  part,  l.e.,  T ■ S(tj,  i),  1 i 1 < y)* 
Also,  let 

Y 

M - u S(t.,  *), 

1-1  1 

with  closure  M and  coaplensnt  Mc  (with  respect  to  [0,  1]). 

Since 

sup  (A  (t;  K)  + n*lt  - «’l(t)]>  I B* 

Osttsl  “ 


It  suffices  to  show  that,  for  any  « > 0, 
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(3.3)  P(sup  (A  (t;  K)  + n*[t  - GT_1(t)])  - B+  > c)  - o(l). 


Ostsl  “ 

We  deal  with  the  eup  In  (3.3)  In  two  parts  eup  _ end  a up  . First, 

t«H  t«MC 

P(eup  {A  (t;  K)  ■+  n*[t  - oT^Xt)}  - E*  > «) 

UM  * n 

i { P(sup  _ (A  (t;  K)  - A (t4J  K)  + n*(t  - <X_1(t)  - 6+J}  > e) 

1-1  t€S(tltl)  11  ° 1 

* I P(«up  _ iw  (tj  K)  - W (t«j  K)l  > e) 

1-1  tcs^.l)  n n 1 

+ { P(*«P  _ |K'(K_1(t))  - K,(K"1(t.))|*|d„  - nJlB-  > c) 

1-1  t«S(t1(t)  1 n Tt  K 

+ { P(eup  _ |K’(k"1(t))K‘1(t)  - K’(K"1(t.))K_1(t.)||B  - 8-1 18„  > c), 

1-1  tcS(t1,Z)  1 1 n K K 

By  an  argument  Identical  to  that  of  Raghavacharl  (1973)  for  the  uniform 
distribution.  It  followa  that  the  final  aummatlcn  In  the  last  Inequality 


may  be  made  arbitrarily  amall  for  appropriate  choice  of  t.  We  now  derive 
alallar  reeulta  for  the  other  two  aummatlona. 

Flrat  conalder  a(t)  - K’ (K_1(t))K"1(t) , 0 s t s 1.  Recalling 
(3.1)  and  (3.2),  define 

(1+0/6.  , 

Z(F.  1)  - (ts  0 « t S 1,  |a(y)  - s(t) I < € 1 X) 


l±  - (2(y,  t) s y c S(tt,  £)),  1 & 1 & y. 


By  assumption  B,  aup  la(t)|  <c  for  aone  c > 0.  Also,  since  GK" 

OStSl 

la  nondecreasing,  S(t^,  l)  Is  an  Interval  and  S(t^,  t)  la  compact. 


mt 
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Further,  2±  is  an  opsn  covsr  for  S(tlf  O.lSiSy.  Therefore,  there 


exists  a finite  subcover  of  Z4  with  cardinality  y i f[2c 4e 
say  {Z(tjj,  4),  1 £ J £ y^>.  Then , for  n sufficiently  large, 

L 

P(“  l(^n  " V/BkI  ’ aup  - l«<t)  - *(t.)|  > e) 

K K t€S(tl,^)  1 

s l P(n*|(0  - 8„)/B  |aup  |s(t)  - *(t.)|  > e) 

J-l  n K K t€Z(t±i>l)  1 

I1  * , (1+0/6,  i 

s l P(n*|(fi  - 0/#J  € 2 2.£"1  > e/4) 

j-1  n a r 


-d+«2)/62i 


5 l P(“  I (B  ~ 

j-1  n 


> 4/ 4« 


(1/Y,) 


;i  (i/6,)  -i  d-w,) 

* c,  l [4«  2 -l  ] 2 


-(1+6  )/6  (1/6  ) (1+0 
* c,  2[2c4e  2 2 + lj[4c  2 4 *)  2 


- o(4  2). 

Note  that  the  fourth  inequality  above  follows  from  Assumption  B(ii). 

*6, 

Since  y - 0(£  ),  for  n sufficiently  large,  the  sum 


l P(»*l <S  - O/BJ 


K)/0Klsup  _ |*(t)  - *(t.)|  > «) 


t«S(tt,4) 


~»°2 

Is  s Y • 0(4  ) and  can  be  made  arbitrarily  small  for  appropriate  choice 

of  4.  A similar  argumsnt  with  *(c)  - K,(K_1(t)),  0 £ t £ 1,  completes  the 


proof  that 


(3.4)  P(aup  _{A  (t;  K)  + n*(t  - GK-1(t) ]}  - E+  > e)  - o(l). 
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Next  we  deal  with  the  eup  part.  Since  GK~^  Is  continuous,  8+  c M 

teMc 

for  every  choice  of  t i.  1 and  sup  [t  - GK  *(t)]  Is  bounded  above  by  a 

teMC 

number  p with  0 < p < 8*.  Choose  n < 8+  - p.  Then 


sup  {A  (t;  K)  + n*(t  - «“l(t>]>-E* 
teMc  n n 

S sup  (A  (t;  K)  - sup  . A (t)  + n^(p  - 8+) 

t«MC  n te6+  n 

s 2 sup  |A  (t;  K)|  + n^(p  - 8+). 

Ostsi  n 


Since  sup  |A  (t;  K) I - 0 (1)  and  n^(p  - 0+)  < -n^n, 

OStsl  n P 

we  have 

(3.5)  P(sup  {d  (t;  K)  + n*(t  - GK_1(t)]>  - E*  > e)  - o(l). 
teM6  n n 

Thus  (3.3)  follows,  completing  the  proof  of  Theorem  1.  □ 

The  proof  of  Theorem  2 Is  now  given.  The  method  also  provides  a 
simpler  proof  for  Theorem  2 of  Raghavachari  (1973). 

PROOF  OF  THEOREM  2.  Write  D ■ max{D  , D+>.  If  S,  and  6,  are  nonempty, 

n n n i z. 

then  6+  ■ 8 • 0~,  0^  ■ 8+,  and  0 2 • 0*".  Therefore,  for  every  a > 0, 

P(D  - n^8  S a)  ■ P(max(D+  - n^6+,  ~(D~  - n^  0)  } £ o ) 
n n n 

P(max(sup  . V°(t)  , -Inf  _V°(t)}  £ a), 
te0  te0~ 

as  n •. 

How  suppose  0^  Is  empty.  Then  8+  < 8.  Since  ©2  cannot  simultaneously 
be  empty,  0”  8 and  ©2  • 0~.  Therefore, 
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Dn  - n+0  - mex{(D+  - nV)  4 n\e+  - 6),  - (o’  - „V)). 

n 

81nC*  Dn  ~ n*e*  “ °p<1>  «nd  6+  - 0 < 0, 

P(D  - n*6  so)-*-  P(inf  V°(t)  s o)  , n > •. 

t«e" 

But  for  eopty,  by  convention, 

aaxUup  V°(t),  -inf  V°(t)}  - inf  V°(t). 
t‘e!  tee2  t<6 

For  e2  oopty,  o similar  argument  givaa  tha  desired  result.  0 

We  now  prove  the  basic  leone  which  wee  used  at  several  points  in 
the  proof  of  Theorem  1.  Por  any  continuous  function  f on  10,  1)  end 
subset  A of  (0,  1J,  end  any  c > 0,  define 

S’(x,  t)  - {y:  0 s t S 1,  |f(x)  - f(y)|  < «} 

end 

S’  - {S(x,  c):  x c A). 

LBH*A.  1st  K be  a eubeet  of  [0,  1J  and  a continuous  flotation  on 
10,  1],  Then  there  exists  a finite  open  cover  of  A,  T(  c s'),  with  car- 
dinality s 2[ [2m/c] ],  where  |f|  s M. 

PROOF.  Let 

D4  - (xi%l<  < f(x)  s *(i  4 l)*).  Oils  [ [2M/« ] ] 
and 

Di4[(2M/<]]  " {xi  “♦<1  ♦ *>«  * *<*)  * -%!<},  Isis  [ [ 2M/« ) ] . 


■ — 
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f- 


i 


Without  lots  of  generality,  It  uy  b«  uiuMd  that  n A * 4, 

Oils  2([2M/«]].  Now  choose  c n A,  0 S 1 s 2([2M/«]],  and  con- 
alder  S'(x.,  c) , 0 i 1 s 2 [ [2M/< ] ] . Suppoae  y e D^.  Then , If  f(x^)  s 0, 

%lc  < f(x)  s (1  + 1)«  and  -*(i  + 1)«  < -f(y)  s -*lc. 

Therefore,  I f (x^)  - f(y)|  & H*.  and  c S'Cx^,  e).  Similarly,  thla  latter 
reault  alao  holda  If  f<Xj)  & 0.  Thua 


Vls2[[2M/«]]Dl  C Ulsis2  [ (2M/<  ] JS ' (*1 ' °* 


But  A c uosis2 ( [2M/e ] ]^1*  Th*r*for* » *>y  construction,  we  have  demonatrated 
the  existence  of  a finite  aubcover  with  the  required  properties.  □ 


4.  Testing  normal  versus  Cauchy.  Suppose  we  want  to  teat 

Hq i F(x)  ■ 4((x  - p)/o),  -•  < x < •»,  for  som  p and  o > 0, 

against  the  fixed  alternative 

H^:  F(x)  ■ C((x  -a)/B),  - «•<  x < *»,  for  som  a end  0 > 0, 

where  C(*)  Is  the  standard  Cauchy  cdf.  In  this  case  p and  o will  be 
estimated  by  the  sample  median  and  a linear  function  of  the  sample  inter- 
quantile range,  respectively.  In  particular,  if  we  define,  for  any 

0 « p < 1,  5 - a + 0C-1(p)  and  t - F_1(p),  where  F is  the  sample  cdf, 

p p n n 

then  (p,  o)  will  be  estimated  by  (p  , 0 ) , where  £ • £,  and 

n n n i 

3n  “ (*3/4  “ ^1/^)/  1.349.  (Note  that  the  Interquartile  range  of  the 
standard  normal  distribution  la  ♦~1(3/4)  - •”l(l/4)  - 1.349.) 
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Since  the  Cauchy  distribution  satisfies  Assumptions  A with  C'(x)  > 0, 
-•  < x < •,  we  havs,  from  Bahadur  (1966),  with  probability  1, 

(4.1)  n*(£p  - Cp)  - n*[p  - Fn<V\{/,(V  + 0(n~*  log  n) . 

This  implies  consistency  of  (0r,  dR) , which  shows  that  ■ a end 

eK  - »<«3/4  - «1,4)]/»(1.349]  ■ B/.6745,  since  the  semi-interquartile 
range  is  the  scale  parameter  of  the  Cauchy.  Therefore, 

K(x)  ■ C(x/.6745),  -•  < x < «. 

Also,  B(li)  follows  from  (4.1)  and  Chebyshev's  inequality  and  B(l)  follows 
from  (4.1)  and  the  Cramdr-Wold  theorem. 

Next  we  need  to  specify  6^,  G^,  0+  and  G . Since  the  minimum  and 
maximum  of  [t  - MC  *(t)]  occur  at  points  -t  ■ K(y) , where  K'(y)  ■ ♦,(y), 
y * 0,  i.e.,  at  y ■ ± .33927,  we  obtain  G^  “ 0 * {.60407}  and 

ex  - g+  - (.39593),  andee^-ee+  - -e”  - 0.0182. 


Combining  this  with  the  fact  that  for  0 s s si,  the  vector 

% 

n ^Wu^8’  *®n  " °K^eK*  ^®n  ~ converges  in  distribution  to 

normal  with  mean  vector  0 and  covariance  matrix 


> 
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we  obtain  that  the  vector  (aup  V°(t) , eup  V°(t))  la  bivariate  normal 

trOj  tce2 

with  mean  0 and  covariance  matrix 


Thua  D+  and  D~,  aultably  normalised,  are  asymptotically  normally  distri- 
buted, while  Dn,  suitably  normalised,  corresponds  to  the  supremum  of 
two  bivariate  normal  r.v.'a  with  sero  mean. 
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